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We report the first observation of tfie decay B" tt^tt", using a 253 fb"'^ data sample coliected at 
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including systematic uncertainties. We also make the first measurement of the direct CP violating 
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Measurements of the mixing-induced CP violation pa- 
rameter sin2(/)i at B factories are in good agreement 
with the Kobayashi-Maskawa (KM) mechanism To 
confirm this theory, one now has to measure the other 
two angles of the unitarity triangle, ^ and (^3 . One tech- 
nique for measuring (j)2 is to study time dependent 
CP asymmetries in i?*^ —> tt+tt" decay, where we have 
recently reported 16] the observation of CP violation and 
evidence for direct CP violation. The extraction of 02, 
however, is complicated by the presence of both tree and 
penguin amplitudes, each with different weak phases. An 
isospin analysis of the tttt system is necessary |7| , and one 
essential ingredient is the branching fraction for the de- 
cay B'^ ^11° 11°. 

QCD-based factorization predictions for B{B^ — > 
7r°7r°) are typically around or below 1 x 10~^ but 
phenomenological models incorporating large rescatter- 
ing effects can accommodate larger values |9|. Evidence 
for B° tt^tt" emerged at the B factories a year 

ago, with a combined value of (1.9 ± 0.5) x 10~^ for the 
branching fraction If such a high value persists, an 



isospin analysis for (j)2 extraction would become feasible 
in the near future. To complete the program, one would 
need to measure both the B° and B° decay rates, i.e. 
direct CP violation. 

In this paper we report the first observation of the de- 
cay B° — > 7r°7r°. We also make a first measurement of the 
direct CP violating asymmetry in this mode. The results 
are based on a 253 fb~^ (275 M BB pairs) dataset col- 
lected with the Belle detector at the KEKB e+e~ asym- 
metric collider ■ KEKB operates at a center-of-mass 
(CM) energy of -/s = 10.58 GeV, corresponding to the 
mass of the T(4S) resonance. Throughout this paper, 
neutral and charged B mesons are assumed to be pro- 
duced in equal amounts at the T(4S), and the inclusion 
of charge conjugate modes is implied, unless otherwise 
specified. 

The Belle detector is a large- solid- angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
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(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside of the coil is instrumented 
to detect mesons and to identify muons (KLM). The 
detector is described in detail elsewhere . Two differ- 
ent inner detector configurations were used. For the first 
sample of 152 million BB pairs (Set I), a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; 
for the latter 123 million BB pairs (Set II), a 1.5 cm ra- 
dius beampipe, a 4-layer silicon detector and a small-cell 
inner drift chamber were used [l^ . 

Pairs of photons with invariant masses in the range 
115 MeV/c^ < m^^ < 152 MeV/c^ are used to form tt^ 
mesons; this corresponds to a window of ±2.5(t about 
the nominal 7r° mass, where a denotes the experimental 
resolution, approximately 8 MeV/c^. The measured en- 
ergy of each photon in the laboratory frame is required 
to be greater than 50 MeV in the barrel region, de- 
fined as 32° < 6*^ < 129°, and greater than 100 MeV 
in the end-cap regions, defined as 17° < 9~f < 32° and 
129° < Oj < 150°, where 9^ denotes the polar angle of 
the photon with respect to the positron beam line. To 
further reduce the combinatorial background, 7r° candi- 
dates with small decay angles (cos 9* > 0.95) are rejected, 
where 9* is the angle between the tt^ boost direction from 
the laboratory frame and one of its 7 daughters in the tt*^ 
rest frame. 

Signal B candidates are formed from pairs of tt*^ 
mesons and are identified by their beam energy con- 
strained mass Mbc = a/ E^lam ~ P*B ^"^^ energy differ- 
ence ls.E = Eg — i^boamj where £^boam denotes the beam 
energy and p*g and E^ are the momentum and energy, 
respectively, of the reconstructed B meson, all evaluated 
in the e+e" CM frame. We require Mbc > 5.2 GeV/c^ 
and -0.3 GeV < A£; < 0.5 GeV. The signal efficiency is 
estimated using GEANT-based [THI Monte Carlo (MC) 
simulations. The resolution for signal is approximately 
3.6 MeV/c^ in A/bc- The distribution in /S.E is asymmet- 
ric due to energy leakage from the CsI(Tl) crystals. If it 
is parameterized by a bifurcated Gaussian, the upper and 
lower resolutions are 46 MeV and 122 MeV, respectively. 

We consider background from other B decays and from 
e+e" {1 = u, d, s, c) continuum processes. A 

large generic MC sample shows that backgrounds from 
h c decays are negligible. Among charmless B decays, 
the only significant background is B^ — > p^7r° with a 
missing low momentum tt^. This background populates 
the negative A£' region, and is taken into account in the 
signal extraction described below. 

The dominant background is due to continuum pro- 
cesses. We use event topology to discriminate signal 
events from this qq background, and follow the con- 
tinuum rejection technique from our previous publica- 
tion j^. We use modified Fox- Wolfram moments 0| 
where the particles in the signal B candidate (category 



s) and those in the rest of the event (category o) are 
treated separately; we also use the missing momentum of 
the event as a third category (category m). Some addi- 
tional discrimination is achieved by considering charged 
and neutral particles in the o category independently, 
and by taking the correlations of charges into account. 
We combine 16 modified moments with the scalar sum of 
the transverse momentum into a Fisher discriminant [l7| , 
and tune the coefficients to optimize the separation be- 
tween signal and background. 

The angle of the i3-meson fiight direction with respect 
to the beam axis {9b) provides additional discrimination. 
A likefihood ratio TZg = Cs/{^s + ^qq) is used as the 
discrimination variable, where denotes the product of 
the individual Fisher and 9b likelihoods for the signal 
and Cqq is that for the qq background. The likelihood 
functions are derived from MC for the signal and from 
events in the Mbc sideband region (5.20 GeV/c^ < Afbc < 
5.26 GeV/c^) for the qq background. 

Additional discrimination between signal and back- 
ground can be achieved by using the Belle standard algo- 
rithm for &-fiavor tagging P, |^ , which is also needed for 
the direct CP violation measurement. The flavor tagging 
procedure yields two outputs: q = ±1, indicating the fla- 
vor of the other B in the event, and r, which takes values 
between and 1 and is a measure of the confidence that 
the q determination is correct. Events with a high value 
of r are considered well-tagged and are therefore unlikely 
to have originated from continuum processes. For exam- 
ple, an event that contains a high momentum lepton (r 
close to unity) is more likely to be a BB event so a looser 
TZs requirement can be applied. We find that there is no 
strong correlation between r and any of the topological 
variables used above to separate signal from continuum. 

We divide the data into r > 0.5 and r < 0.5 bins. The 
continuum background is reduced by applying a selection 
requirement on TZ^ for events in each r region of Set I and 
Set II according to t he figure of m erit (FOM). The FOM 
is defined as iV^'^P / ^iVj^Tlvff , where iV,'=''P and N";^^ 
denote the expected signal, assuming the branching frac- 
tion B = 2 X 10~^, and background yields obtained from 
MC and sideband data, respectively. A typical require- 
ment suppresses 97% of the continuum background while 
retaining 53% of the signal. 

The signal yields are extracted by applying unbinned 
two-dimensional maximum likelihood (ML) fits to the 
(Mbc, Ai?) distributions of the B and B samples. The 
likelihood is defined as 

C = CXp (- J2 ^-fc.j) Ns,k,JVs,k,J,^) (1) 

s,k,j i s^k^j 

where 

Vs,k,J,^ = P.^fcj(Mbc»,Ai;,), (2) 

and s indicates Set I or Set II, k distinguishes events in 
the r < 0.5 or r > 0.5 bins, i is the identifier of the i-th 
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FIG. 1; Result of the fit described in tiie text. (Left) Mhc pro- 
jection for events that satisfy -0.2 GeV < AB < 0.05 GeV; 
(right) AS projection for events that satisfy 5.27 GeV /(? < 
Afbc < 5.29 GeV/c^. The solid lines indicate the sum 
of all components, and the dashed, dotted and dot-dashed 
lines represent the contributions from signal, continuum, and 
P^tt", respectively. 

event, Ps.k.jiM^c, ^E) are the two-dimensional proba- 
bility density functions (PDFs) in Mbc and AE for the 
signal and background components, Nj is the number of 
events for the category j, which corresponds to either 
signal, qq continuum, or background from —> p^n^ 
decay. 

The PDFs for the signal and for 5+ — > p+7r° are 
taken from smoothed two-dimensional histograms ob- 
tained from large MC samples. For the signal PDF, 
discrepancies between the peak positions and resolutions 
in data and MC are calibrated using — > tt^tt^ and 
B+ i^+7r-7r°)7r+ decays. The difference is 

caused by imperfect simulation of the 7r° energy reso- 
lution while the effect of the opening angle distributions 
can be neglected. The invariant mass distribution for 
the is fitted with an empirical function for data and 
MC, and the observed discrepancies in the peak position 
and width are converted to the differences in the peak 
position and resolution for AE in the signal PDF. We 
require the decay products to lie in the same momen- 
tum range as the tt'^s from B 7r°7r°. To obtain the 
two-dimensional PDF for the continuum background, we 
multiply the PDF for AE, which is modeled with a lin- 
ear function, with the PDF for Mbc, for which we use 
the ARGUS function In the fit, the shapes of the 

signal and B^ p+Tr^ PDFs are fixed, with the normal- 
ization for B^ p+7r° fioated; all other fit parameters 
are allowed to float. The flt results are shown in Fig.^ 

The obtained signal yield is 81.8l]^g g with a statisti- 
cal significance (S) of 6.1, where S is defined as S — 
yJ—2\n{La/LN^), and Co and Cn, denote the maximum 
likelihoods of the fits without and with the signal compo- 
nent, respectively. The relative yields in sets I and II are 
consistent with the expectation based on their relative 
luminosities. We vary each calibration constant for the 
signal PDF by ±lcr and obtain systematic errors from 



the change in the signal yield. Adding these errors in 
quadrature, the significance including systematic uncer- 
tainties is reduced to 5.8tT, which corresponds to the first 
observation of B^ tt^tt". 

In order to obtain the branching fraction, we divide 
the signal yield by the reconstruction efficiency, measured 
from MC to be 12.9%, and by the number of BB pairs. 
We consider systematic errors in the reconstruction effi- 
ciency due to possible differences between data and MC. 
A 4.2% systematic error is assigned for the uncertainty 
in the efficiency for the track multiplicity requirement. 
This is determined by varying the multiplicity distribu- 
tion of signal MC. We assign a total error of 6% due to 
TT^ reconstruction efficiency, measured by comparing the 
ratio of the yields of the t] Tr^Tr^Tr" and ry — > 77 de- 
cays. The experimental errors on the branching fractions 
for these decays are included in this value. We check 
the effect of the continuum suppression using a control 
sample of 5+ £'''(—> if +7r~7r°)7r+ decays; the TZg re- 
quirements has a similar efficiency for the MC control 
sample and for signal MC. Comparing the TZg require- 
ment on the control sample in data and MC, a systematic 
error of 1.8% is assigned. We check for a possible pile-up 
background due to hadronic continuum events that con- 
tain energy deposits from earlier QED interactions. Such 
a background may peak in Mbc, however, the showers 
from the QED interaction can be identified from timing 
information recorded in the ECL. For Set II, it is pos- 
sible to remove these events using this information and 
determine the change in event yield. We conservatively 
estimate a systematic uncertainty of 10.3% for this off- 
time QED background. Finally, we assign a systematic 
error of 1.1% due to the uncertainty in the number of 
B'B pairs (274.8 ± 3.1) x 10*^, and obtain a branching 
fraction of 

B{B' - tt^tt") = (2.31°:™) X 10-^. 

The result is stable under variations of the TZg cut. 

Having observed a significant signal, we utilize the 
B^ / B^ separation provided by the flavor tagging to mea- 
sure the CP asymmetry. Equation (2) is replaced by 

Vs,k,j = ]^[\-q^■Acp'LJ]Ps±^J{Mhc^.AEi), (3) 

where q indicates the B meson flavor, B{q = -1-1) or 
B{q = —1), Acp'i,j is the effective charge asymme- 
try, where Acp'i,j = AcPj{l - 2x<i)(l - 2w;). Here 
Xd = 0.186 ± 0.004 is the time-integrated mixing 
parameter and wi is the wrong-tag fraction. For qq con- 
tinuum, Xd and wi are set to zero. The tt^'tt" sample is 
divided into six r-bins, and the r-dependent wrong-tag 
fractions, wi {I = 1, . . . , 6), are determined using a high 
statistics sample of self-tagged B^ D^*'>^n^ , D*^ 
and D*~£^v events and their charge conjugates j20| . The 
total number of signal events is fixed to the yield obtained 
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FIG. 2: Mbc and AiJ distributions with projections of the 
fit superimposed. The distributions are shown separately for 
events tagged as B° (left) and (right). 



branching fraction significantly higher than factorization 
predictions. We obtain 81.8^i[gg signal events with a sig- 
nificance of 5.8 standard deviations (a) including system- 
atic uncertainties. The branching fraction is measured 
to be (2.3t°:^t° J) X 10-*'. This result is consistent with, 
and supersedes, our previous result. We have also made a 
first measurement of the direct C P violating asymmetry. 
The large branching fraction for tt^tt", together 

with the measurements of its direct CP violating asym- 
metry Acp, will allow a model- independent extraction 
of the CKM angle (/)2 from measurements of the B ^ tttt 
system in the near future. 
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from the branching fraction measurement. The relative 
fractions of signal events, qq and p*7r° background events 
in the different r bins are also fixed. 
Defining the direct CP asymmetry as 



Acp = 



_ N{B ^ /) - NiB ^ /) 



NiB^ f)+NiB^ fy 



(4) 



the result is Acp = OAAZo%2 ± O-l'^- Systematic errors 
are estimated by varying the fitting parameters by zblcr. 
Including the result of a null asymmetry check with the 
same analysis procedure for the B D{KTTn^)TT control 
sample, the total systematic error is ±0.17. To illustrate 
this asymmetry, we show the results separately for B'^ 
and B^ tags in Fig. [3 While not significant, the method 
already gives constraints on HJ] . 

Our results confirm the previous evidence [lol ll]| 
and establish the decay 5° tt^tt^. Since the ob- 
served branching fraction is much larger than predic- 
tions based on QCD factorization recent theoreti- 
cal discussions have focused on the possibility of an en- 
hanced color-suppressed amplitude, together with a siz- 
able strong phase Other color-suppressed modes 
such as B° £)°7r° and B^ — > p°7r° have also been mea- 
sured |U 0| at rates considerably higher than factoriza- 
tion predictions 0, . In addition, the recent evidence 
for large direct CP violation in B^ tt+tt" |^ and 
K+TT' modes '2^ disagrees with QCD based fac- 
torization predictions. Some effect beyond factorization 
appears to be present in charmless two-body B decays. 

In conclusion, we have observed the B^ tt^tt" decay 
mode in a data sample of 275 million BB pairs with a 
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